arXiv:1504.05954v2 [cond-mat.mtrl-sci] 13 Sep 2015 


Effect of inhomogeneities and substrate on the dynamics of the metal-insulator 

transition in VO2 thin films 

M. Rodriguez-Vega/ M. T. Simons/ E. Radue/ S. Kittiwatanakul/ 

J. Lu/ S. A. Wolf/’^ R. A. Lukaszew/ I. Novikova/ and E. Rossi^ 

^Department of Physics College of William and Mary Williamsburg VA, 23187 
^Department of Material Science University of Virginia Charlottesville VA, 22904 
^Department of Physics University of Virginia Charlottesville VA, 22904 
(Dated: September 15 , 2015 ) 

We study the thermal relaxation dynamics of VO2 films after the ultrafast photo-induced metal- 
insulator transition for two VO2 him samples grown on AI2O3 and Ti02 substrates. We hnd two 
orders of magnitude difference in the recovery time (a few ns for the VO2/AI2O3 sample vs. hundreds 
of ns for theV02/Ti02 sample). We present a theoretical model to take into account the effect of 
inhomogeneities in the hlms on the relaxation dynamics. We obtain quantitative results that show 
how the microstructure of the VO2 him and the thermal conductivity of the interface between the 
VO2 him and the substrate affect long time-scale recovery dynamics. We also obtain a simple 
analytic relationship between the recovery time-scale and the hlm’s parameters. 


I. INTRODUCTION 

Vanadium dioxide (VO2) undergoes a metal-insulator 
transition (MIT) around room temperature^ enabling 
a wide range of potential applications. It has re¬ 
cently been shown that it is possible to photo-induce 
the insulator-to-metal transition in VO2 in the sub¬ 
picosecond timescale^^®. This finding makes VO 2 a ma¬ 
terial of great interest for electronic and photonic appli¬ 
cations, such as ultra-fast switches or transistors. The 
realization of V02-based switches requires the ability to 
control the VO2 MIT dynamics using external fields, as 
well as a better understanding of the recovery mecha¬ 
nisms after the external field is turned off and the ma¬ 
terial returns to its normal state. The mechanism by 
which the photo-induced insulator-to-metal takes place 
in VO2 is still not fully understood due to the complex¬ 
ity of the electronic behavior of VO2 arising from the 
presence of strong electron-lattice coupling and electron- 
electron interactions®^^^. As a result, VO2 is a unique 
material of great fundamental and practical interest. 

At low temperatures (T < 340 K) the VO2 lattice 
has a monoclinic structure, whereas at high temperatures 
(T > 340 K) it has a tetragonal structure. This differ¬ 
ence in lattice structure is reflected in the band struc¬ 
ture: VO2 is an insulator in the monoclinic phase and 
a metal in the tetragonal phase. This simple picture is 
complicated by the fact that in VO2 electron-electron cor¬ 
relations are very strong and can provide an important 
contribution to the localization of the electronic states 
via the Mott mechanism^°di,i3,i4^ appears that a full 
account of the MIT must take into account the inter¬ 
play of the lattice dynamics and the electron dynamics 
driven by strong electron-electron interactions. This is 
a fascinating and extremely challenging problem that in 
addition is complicated by the unavoidable presence of 
inhomogeneities 1^46^ 

Several works®’^^^^® have investigated the short 
timescale dynamics after the photo-induced transition. 


In particular. Ref. 26 presented a comparison of the 
long timescale recovery dynamics between VO2 films on 
a crystal substrate or a glass substrate and found that 
the recovery time for the films on the glass substrate was 
much longer than for the films on a crystal substrate. 
The recovery time was modeled using the heat equation 
to describe the heat flow across the interface between the 
VO2 film and the substrate. The difference in the char¬ 
acteristic time between the two types of substrates was 
attributed to the fact that the thermal conductivity of 
the interface was expected to be much smaller for glass 
substrates than for crystal substrates. 

In this work we present a theory to properly take into 
account the effect of inhomogeneities on the recovery dy¬ 
namics of VO2 films. Our theory describes simultane¬ 
ously: (i) the profile of the reflectivity across a thermal 
induced MIT; (ii) the long timescale recovery dynamics of 
the reflectivity after a photo-induced insulator-to-metal 
transition; (iii) the observed difference of two orders of 
magnitude between samples with different substrates. In¬ 
homogeneities are due to the fact that the film is com¬ 
prised of grains with different sizes and different local 
properties, such as strain^^’^s and chemical composition. 

The presence of inhomogeneities induces a distribu¬ 
tion of values for the transition temperature Tc within 
the film. To take this into account we derive a gener¬ 
alized heat equation that includes the fact that during 
the recovery from the photo-induced insulator-to-metal 
transition, at any given time a fraction of the sample 
is undergoing the metal-to-insulator transition, another 
fraction is still cooling in the metallic phase, and an¬ 
other fraction is already cooling in the insulating phase. 
A key ingredient of the generalized heat equation is the 
correct description of the time evolution of the fraction 
of the sample that is metallic, insulating, or undergo¬ 
ing the phase transition. We then use our theoretical 
model to obtain the scaling relation between the charac¬ 
teristic recovery time r and the parameters of the films. 
Our theoretical model, and the underlying assumptions. 
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are strongly supported by our experimental results. Dif¬ 
ferently than in Ref. 26 our VO 2 films have all crystal 
substrates, no glass. Yet, we Hnd that r can be more 
than two orders of magnitude different depending on the 
crystal substrate, Ti 02 or AI 2 O 3 , Fig. 1 . 

The generalized heat equation, Eqs. (9), which prop¬ 
erly takes into account the effect of the films inhomo¬ 
geneities on the recovery dynamics is the main result of 
our work. Our theory allows the description of the re¬ 
covery dynamics consistently with the measurements ob¬ 
tained for the thermally driven MIT. The scaling between 
the characteristic recovery time t and the parameters of 
the film is another important result of our work. 




Time (ns) 

FIG. 1. Relative change in reflectivity (AR/R) for the VO 2 
film on (a) AI 2 O 3 substrate and (b) Ti 02 substrate as a func¬ 
tion of time after the MIT is induced at time t = 0 by a 
strong ultrafast pump pulse. The values of the pump flu- 
ence are shown in the legend, and the sample temperature 
is set to 311 K in (a) and 280 K in (b), which correspond 
to approximately 30 K below the critical temperature for 
thermally-induced MIT for each sample. 

Our work is relevant to the more general problem of 
how spatial inhomogeneities affect a first order phase 
transition. The ability of our treatment to contribute 
to this general problem relies on the fact that in VO 2 the 
two phases across the first-order phase transition have 
very different electronic properties (metallic vs. insulat¬ 
ing behavior) that allows us to get an accurate phase 


mapping, via optical reflectivity measurement, of the 
time evolution of the metallic (insulating) fraction and, 
indirectly, of the spatial inhomogeneities present during 
the transition. 

The work is organized as follows. Sec. II describes 
the experimental arrangements to measure the optical 
reflectivity time-evolution. The details of the theoreti¬ 
cal model that we use to characterize the distribution of 
the films’ inhomogeneities and the long-time dynamics of 
the reflectivity after a photo-induced insulator-to-metal 
transition are presented in Secs. Ill and IV, respectively. 
In Sec. V we demonstrate how the variations in statistical 
properties of the two films result in a significant difference 
in the relaxation timescales, and in Sec. VI we provide 
our conclusions. 


II. EXPERIMENTAL SETUP 

In our experiments we studied two VO 2 thin-film sam¬ 
ples, both of which were produced using reactive-bias tar¬ 
get ion beam deposition (RBTIBD)^®. One sample was 
grown on 0.5 mm thick C-AI 2 O 3 , and the thickness of the 
VO 2 film was 80 nm. The other sample was grown on 
a 0.5 mm thick Ti 02 (011) substrate, and was measured 
to be 110 nm thick. X-ray diffraction (XRD) evaluation 
of both films showed them to be crystalline, and detailed 
characterization information is available in previous re- 
ports^®’^^. 



VO2 Sample 


FIG. 2. Schematic of the ultrafast pump-probe setup. BS is 
an 80/20 beam splitter. 

For the photo-induced insulator-to-metal transition ex¬ 
periments we used an ultrafast laser system (Coherent 
Mantis oscillator and Legend Elite regenerative ampli¬ 
fier) with approximately 100 fs pulses with a central 
wavelength at 800 nm and a repetition rate of 1 kHz. 
The properly attenuated output of the laser was split 
into strong pump pulses and weaker probe pulses using 
a beam splitter in a standard pump-probe configuration, 
shown in Eig. 2. The more powerful pump beam, focused 
to a 180 /im diameter spot on the surface of the sample, 
was used to induce the insulator-to-metal transition, and 
its fluence was controlled using a variable neutral-density 
filter (VF). The fluence of the probe beam was further 
attenuated to a value well below the insulator-to-metal 
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threshold ((^probe < 100 /rJ/cm^), and we used its re¬ 
flectivity from the sample to monitor the instantaneous 
optical properties of the VO2 him. The probe pulses were 
directed along a variable delay stage to accurately control 
the relative timing between the pump and probe pulses 
by up to 4 ns with a few fs precision. The probe beam 
was focused on the sample at the same spot as the pump 
beam, using a shorter focal length lens. When tuned 
to the center of the pump beam focal spot, the smaller 
probe beam diameter (90 /rm) ensured probing a region 
of uniform pump intensity. 

The reflected probe power was measured using a sili¬ 
con photodetector, and further analyzed using a lock-in 
amplifier. To minimize the effects of probe pulse insta¬ 
bilities, as well as long-terms drifts due to environmental 
changes, we report the relative change in probe reflection 
AR/R with the pump beam on or off. 


Variable ND filter 



VO2 Sample 

FIG. 3. Schematic of the experimental setup using a 
continuous-wave probe laser. 

Notably the MIT relaxation of the V02/Ti02 sample 
was not measurable with the femtosecond probe, as its 
characteristic decay time exceeded the 4 ns maximum 
pulse separation, determined by the length of the de¬ 
lay stage. To measure the relaxation of the metallic VO2 
grown on the rutile sample we modified our experimental 
setup by replacing the femtosecond probe pulses with a 
continuous-wave (CW) diode laser operating at 785 nm 
and a fast photodiode (measured response time of ap¬ 
proximately 10 ns), as shown in Fig. 3 . This detection 
method allowed us to measure changes in reflectivity for 
times longer than « 20 ns after the insulator-to-metal 
transition, that were inaccessible with the femtosecond 
probe arrangement. 

Figure 1 shows sample measurements of both the 
VO2/AI2O3 and V02/Ti02 films, using the femtosecond 
and CW probe arrangements respectively. The overall re¬ 
flectivity depends on the refractive index of both the film 
and the substrate, and the refractive indices of Ti02 and 
AI2O3 are different. Because it is easier to average the 
CW laser reflection signal, the curves for V02/Ti02 are 
smoother than the curves for the VO2/AI2O3. The rutile 
reflection spectra recorded using the ultrafast probe had 
the same noise as for the sapphire samples, indicating 



FIG. 4. Dependence of metal state decay constant r on the 
laser pump fiuence and substrate temperature. Dots repre¬ 
sent experimental data, and lines corresponds to the results 
of the theoretical calculations. The initial temperature Ta for 
both sample substrates was approximately 30 K below their 
respective MIT critical temperatures. 

that the differences in the noise are due to differences in 
the probes, not in the samples. 

For values of the pump fiuence higher than a thresh¬ 
old, which depends on the substrate temperature, we can 
see that the reflectivity, soon after the pump pulse, re¬ 
mains almost constant for some time, i.e. its dynam¬ 
ics exhibits a “fiat” region, see in particular Fig. 1 (b). 
The observed “flattening” of the curves is due to the 
pump pulse heating the sample to a temperature above 
the threshold value for the thermally-induced insulator- 
to-metal transition^^’^^. In this case the reflectivity stays 
unchanged at the level corresponding to a fully metallic 
phase until a non-negligible fraction of the sample cools 
down to the transition temperature. For all experimental 
curves only the later exponential part of the measured re¬ 
flectivity was included into the fitting thermal relaxation 
time analysis. 

The analysis of the relative reflectivity for both VO2 
samples demonstrate that after the initial rapid change 
during the ultra-fast insulator-to-metal transition, its 
time evolution during the recovery is well fitted by a sin¬ 
gle exponential function with a recovery time constant 
r: 

Rmit)=Ri + iRo-Ri)e-^/^, ( 1 ) 

where Ri corresponds to the reflectivity in the insulating 
phase, and Rq corresponds to the reflectivity at t = 0 s. 
The results of such measurements are shown in Fig. 4 : 
for VO2/AI2O3 films we obtained values of r of the or¬ 
der of few nanoseconds, whereas it took the V02/Ti02 
sample a few hundred nanoseconds to relax back to the 
insulating state. This two orders of magnitude difference 
in the recovery times was even more surprising consider¬ 
ing that the characteristic times for the transition itself 
were quite similar, as demonstrated in previous studies^®. 
In the discussion below we demonstrate that the relax¬ 
ation dynamics strongly depend on the microstructure of 
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the VO 2 films which in turn is strongly influenced by the 
properties of the substrate and their interface. Figure 4 
also reveals that the rate of thermal relaxation for both 
samples increases with higher pump power. 


III. THEORETICAL MODELING OF 
INHOMOGENEITIES 


In order to take into account the effect of the inhomo¬ 
geneities on the MIT dynamics the first step is to char¬ 
acterize them. To do this we can use the profile of the 
reflectivity across the thermally induced MIT. The dot¬ 
ted lines in Figures 5 (a), and (b) show the measured 
reflectivity as a function of temperature across the ther¬ 
mally induced MIT for a VO 2 film grown on sapphire and 
Ti 02 , respectively. The temperature driven MIT in VO 2 
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FIG. 5. Evolution of the reflectivity across the thermally 
induced MIT for the case of sapphire and rutile substrates 
normalized to the average critical transition temperature. 
The open circles (red) correspond to the measured reflectiv¬ 
ity in the heating branch, the solid circles (blue) correspond 
to the measured reflectivity in the cooling branch, and the 
solid curve corresponds to the theoretical result. For rutile 
substrate (Tc) = 314.0 K, and for the sapphire substrate 
(Tc) = 340.1 K. 


is a first-order transition. In the ideal case the reflectiv¬ 
ity is expected to exhibit a finite, step-like, change at the 
critical temperature Tc, at which the sample goes from a 
low-temperature insulating state to a high-temperature 
metallic state. In thin films, however, the optical reflec¬ 
tivity smoothly changes from the value corresponding to 
the insulating phase {Ri) to the value characteristic to 
the metallic phase (Rm) as the temperature increases, as 
Fig. 5 illustrates. For our samples the hysteresis loop is 
very narrow^^. The fact that the MIT takes place over 
a range of temperatures implies that different regions of 
the sample have different values of T^. This is different 
from the case of an ideal, homogeneous, system for which 
the whole sample exhibits the coexistence of metallic and 
insulating regions only for T = T^. As a consequence the 
MIT in the films is characterized not by a single critical 
temperature but by a distribution P{Tc) of critical tem¬ 
peratures. This is due to the fact that the VO 2 films are 
inhomogeneous: they are formed by crystal grains with 
different local properties. Different grains in general have 
different sizes, slightly different stoichiometry, and expe¬ 


rience different local strains. It is very challenging to 
characterize the distribution of all the local properties 
that can affect the transition temperature of each grain. 
However, for our purposes, we only need P{Tc) and, as we 
show below, this can obtained directly from the profiles 
of R{T) without having to characterize the distribution 
of the local properties affecting T^. Let rij be the frac¬ 
tion of the sample in the insulating phase. At a given 
temperature T we have: 

poo 

m{T) = P{Tc)dT,. (2) 

Let ? 7 m(F) = 1 — rjiiT) be the fraction of the film in the 
metallic phase. To obtain the evolution of r]i{T) across 
the MIT, and therefore P{Tc), considering that changes 
in the fraction of the film that is metallic (insulating) are 
the dominant cause of changes in the reflectivity, we can 
use a two-fluid effective medium theory (EMT)^^“^®. In 
the EMT the inhomogeneous system is replaced by an 
effective homogeneous medium having the same, bulk, 
electric properties. Let em, £/ be the dielectric constants 
(at the probing frequency) of VO 2 in the metallic and in¬ 
sulating phase respectively. Then, the dielectric constant 
of the effective medium, cemt-, is given by the following 
equation: 

— (emt) _ — (emt) _ Q 

f-EMT + g{f-I — ^EMt) ^EMT + gi^M — ^EMt) 

( 3 ) 

In Eq. (3) g is a factor that depends on the shape of the 
grain. Without loss of generality we set g = 1/3. Let 
n and k be the real and imaginary parts respectively of 
the index of refraction, so that for the effective medium 
n + ik = y'e emt and therefore 


R = 



cos 00 + 




sin^ 60 


( 4 ) 


where Oq ~ 15° corresponds to the probe incidence angle. 
Given our experimental setup we can reliably obtain the 
imaginary part of the index of refraction by measuring 
the absorption. For this reason we set the value of the 
imaginary part of the complex index of refraction /cjvfj 
{ki) for the metallic and (insulating) phase to the mea¬ 
sured values, consistent with the values reported in the 
literature^®’^^, and then use Eq. (4) and the measured 
value of Rm (Ri) in the metallic (insulating) phase to 
fix the corresponding value of um {ni) (see Table I). 

Using Eqs. (2)-(4) we can obtain the profile of R{T) 
across the MIT for a given P{Tc). Assuming P{Tc) to be 
a Gaussian distribution, by fitting the measured profile 
of R{T) to the one obtained using Eqs. 2-4 we can obtain 
the average value of the critical temperature (Tc) and 
its standard deviation (Tt„. For V 02 /Ti 02 samples we 
find (Tc) = 314 K, = 2.6 K, for VO 2 /AI 2 O 3 samples 
(Tc) = 340 K, ax^ = 8.8 K. The solid lines in Fig. 5 show 
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the profiles of R{T) obtained using Eqs. ( 2 )-( 4 ) and the 
above values for (Tc) and 

The difference in the value of Tc between V02/Ti02 
and VO2/AI2O3 samples can be attributed to the fact 
that Ti02, having a rutile structure, might induce strains 
into the VO2 film that should favor the metallic phase of 
VO2. In general, strain effects are expected to play an 
important role in the physics of the MIT phase transi¬ 
tion of VO2 films. In our approach such effects enter 
indirectly, via the form of the probability distribution 
P{T(.), and the value of the thermal conductivity of the 
interface between the VO2 film and the substrate. 

As we discuss in the following section, for our the¬ 
oretical treatment of the recovery dynamics over long 
timescales of VO2 films the knowledge of P{Tc), i.e., (T^) 
and (TTcj is ah that is needed. As mentioned before the 
fact that (TTc is nonzero is due to inhomogeneities, of dif¬ 
ferent nature, present in the VO2 film. It is practically 
impossible to know the distribution in the films of the 
properties affecting Tc- However, it is interesting to con¬ 
sider the limit in which the grain size D is the dominant 
property affecting T^. The reason is that in this limit it is 
possible to extract, using strong and fundamental argu¬ 
ments, the distribution, P{D), for the grain size. In par¬ 
ticular, it is possible to obtain the average grain size, (D), 
and its standard deviation, quantities that are of great 
practical interest. (D) can be compared to estimates ob¬ 
tained using more direct experimental techniques, such 
as XRD. In the remainder of this section we use the ex¬ 
perimental results for R{T) across the MIT to extract 
{D) and its standard deviation. 

Theoretical and experimental evidence^® indicates that 
for thin films the distribution P{D) of the grain size D 
typically follows a logarithmic-normal distribution. 




In D/i) 
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2a2 
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Using Eqs. ( 2 )-( 7 ), by fitting the measured profile of 
R{T) across the MIT, we can obtain P{D) and therefore 
{D) and its standard deviation. Figures 6 (a), (b) show 
the profiles of P{D) used to obtain the good theoreti¬ 
cal fit to the evolution of R{T) shown in Fig. 5 . These 
profiles correspond to {D) = 64.7 nm ud = 38.5 nm for 
VO2/AI2O3 samples, and {D) = 17.4 nm ax, = l.I nm 
V02/Ti02 samples. It is interesting to compare the val¬ 
ues of {D) obtained using this approach to the values 
estimated using XRD. From XRD data^® we estimated 
{D) « 45 nm for VO2/AI2O3 and {D) « 13 nm for 
V02/Ti02 (see Table I). These values are in remarkable 
semi-quantitative agreement with the values extracted 
from the profiles of R{T) across the MIT suggesting that 
the assumption that the grain size is the dominant prop¬ 
erty affecting the local value of might be qualitatively 
correct. It is therefore interesting to obtain the profiles 
of P(Tc) corresponding to the distributions of grain sizes 
shown in Figs. 6 (a), (b). Such profiles are shown in 
Figs. 6 (c), (d). The evolution of 777 (T) across the MIT 
obtained using these profiles is shown in Fig. 7 . 



VO 2 /Ti02 

VO 2 /AI 2 O 3 

(Tc) 

314.0 K 

340.1 K 

(TT„ 

2.6 K 

8.8K 

(D) 

17.4 nm 

64.7 nm 


13 nm 

45 nm 

(70 

1.1 nm 

38.5 nm 

71m + 

1.53-b jO .8 

1.49-bi0.65 

Rm 

0.14 

0.11 

Hi + ikj 

3.03-ti0.57 

2.60 -I- i0.60 

Ri 

0.28 

0.23 

(7K 

1,100 W/(K cm^) 13, 000 W/(K cm^) 


TABLE I. Comparative table between V 02 /Ti 02 , and 
VO 2 /AI 2 O 3 sample parameters. 


In Eq. ( 5 ) D is the effective diameter of a grain, D is the 
grain size (diameter) such that InZ? = (InD), and a is 
the standard deviation of ln(D). 

From general and fundamental arguments®®^^^ we 
have: 

= P""' (' - dw) • 


where is the bulk transition temperature and Dq, 

equal to 2 nm in our case, is the grain’s diameter be¬ 
low which the grain is so small that is not possible to 
unambiguously identify its crystal structure. We set 
^qbuik) _ ggg ihaA, is the temperature above which the 
VO2/AI2O3 samples is completely metallic. This value is 
higher than the value of bulk VO2 due to the strain ex¬ 
perienced by the films^^’^®. The relation between P{D) 
and P{Tc) is given by: 


p(re) = p(D(T,))—. 


( 7 ) 


Our analysis suggest that the R{T) profiles could be an 
indirect method to characterize the distribution of grain 
sizes in VO2 films, a very challenging quantity to obtain 
using direct imaging experiments. 


IV. THEORETICAL MODELING OF THE 
RELAXATION DYNAMICS OF THE MIT 


In our experiment the VO2 films have a thickness d 
equal to or smaller than 110 nm (see Fig. 8), which is 
comparable with the laser 1/e penetration depth 6 ~ 
110 — 130 nm^®. Thus, we can assume that the pump 
pulse heats the film uniformly throughout its thickness. 
To describe the heat transfer process between the film 
and the substrate, we assume the temperature to be uni¬ 
form throughout the film for all times. Effectively, given 
these conditions, the heat transfer problem becomes a 
one-dimensional problem, and the equation for the rate 
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FIG. 6 . (a) and (b) show the grain size distributions normal¬ 
ized to the average grain size for sapphire {{D) — 64.7 nm) 
and rutile {{D) = 17.4 nm) substrate respectively, (c) and 
(d) show the critical temperature distribution normalized to 
the average critical temperature for sapphire ((Tc) = 340.1 
K) and rutile {{Tc) = 314.0 K) respectively. The bulk critical 
temperature is taken to be _ ggg 


a) b) 



T/{Tc) 


TUT) 


FIG. 7. Evolution of the insulating partial volume rji across 
the thermally induced MIT for case of (a) sapphire and (b) 
rutile substrates. For rutile, {Tc) = 314.0 K, and for sapphire 
{Tc) = 340.1 K. 
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L 

FIG. 8 . Sketch of the heterostructure considered in this work. 
It is composed of a vanadium dioxide (VO 2 ) thin-film de¬ 
posited on top of a substrate. The substrates considered in 
this work are titanium dioxide (Ti 02 ), and aluminum oxide 
(AI 2 O 3 ). For V 02 /Ti 02 d = 110 nm while for VO 2 /AI 2 O 3 
d = 80 nm. For both substrates, L = 0.5 mm. 
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of change of the heat {Q) in the film takes the form: 

=A X d X { piCirii{Tf) + pMCM'r]M{Tf) 

+ L{Tf)P{Tf)pa.) ^ , ( 8 ) 

where Tf is the film temperature, A is the area of the 
film, pi {pm) is the density in the insulating (metallic) 
phase, pav = {pi + Ci (Cm) is the heat capacity 

in the insulating (metallic) phase, L is the specihc heat 
and P{Tf)dTf is the fraction of the sample undergoing 
the MIT in the time interval dt during which the film 
temperature is in the interval [Tf,Tf + dTf], Here P{Tf) 
is the distribution of critical temperatures due to the in¬ 
homogeneities that we have obtained in the previous sec¬ 
tion. Using Eq. ( 2 ) and the fact that tjm = (1 — Pi) we 
can rewrite Eq. (8) in a form that more explicitly shows 
the effect due to the inhomogeneities, i.e. the fact that 
the MIT is not characterized by a single Tc, but by a 
distribution of Tc’s: 


^=Axdx 

dt 


PmCm + {piCi — pmCm) / P{Tc)dTc 


+ L{Tf)P{Tf)pa. 


dJl 

dt 


( 9 ) 


Equation ( 9 ) is the main result of our work: it allows to 
properly take into account the effect of inhomogeneities 
on the long timescale dynamics across a first order phase 
transition. The key quantity entering Eq. ( 9 ) is the dis¬ 
tribution P{Tc) that, as we have shown in the preceding 
section, can be obtained from the prohle of R{T) across 
the thermally activated MIT. Our work is the first to 
combine the information from the thermally activated 
MIT to obtain a physically accurate heat equation to de¬ 
scribe the recovery dynamics after a photo-induced MIT. 

For the latent heat we have^®^"*'^ 


L = 


Tc 

rp(bulk) ■ 
c 


( 10 ) 


where T*^*’'^**') is the value of the specific heat for bulk 
VO2. Given Eq. (6), Eq. ( 10 ) implies L = ~ 

Do/D). 

The rate of change of heat in the him given by Eq. (8) 
must be equal to the heat current {Jq) across the inter¬ 
face between the him and the substrate: 


jQ = -aKA{Tf-Ts{d)) ( 11 ) 


where aK is the Kapitza constant characterizing the ther¬ 
mal conductivity of the interfaces"^^^^® , and Ts{d) is the 
temperature of the substrate at the surface facing the 
VO2 him. Combining Eq. ( 9 ) and Eq. ( 11 ), for Tf we 
obtain the equation: 


PmCm + {piCi — pmCm) / P{Tc)dTc 

•'Tf 


+ L{Tf)P{Tf)pa. 


dTl 

dt 


-^{Tf-T,{d)). 


( 12 ) 
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In Eq. (12) the only undetermined quantity is ax- We 
fix ax by fitting the theoretically obtained time evolu¬ 
tion of R(t) to the one measured experimentally, for fixed 
experimental conditions such as the temperature of the 
substrate and the pump fiuence. The robustness of the 
theory presented is evidenced by the fact that, the same 
fixed value of ax provides a good agreement between the 
theoretical and the experimental results for a broad range 
of experimental conditions. 

To completely define the problem we need to supple¬ 
ment Eq. (12) with proper boundary conditions. The 
temperature distribution within the substrate, Ts(z,t), 
satisfies the diffusion equation: 

dTfzf) ^ ks d'^Tfzf) 

dt C,ps dz^ ^ ’ 

where /Cg, Cg, Ps are the thermal conductivity, heat ca¬ 
pacity, and mass density, respectively, of the substrate. 
The bottom of the substrate, for which z = L (see Eig. 8), 

/ D\ 

is kept at a fixed temperature Tg . At the film/substrate 
interface the heat transferred from the film must be equal 
to the heat current kgdTg/dz\z=d- We then have that the 
boundary conditions completing Eq. (13) are: 

T,(z = L,f) = Ti^); (14) 

= -ax{Tf{t)-Tg{z = d,t)). (15) 

z—d 


dTg{z,t) 

dz 


Equations (10),(12)-(15), combined with knowledge of 
the distribution P{Tc) completely define the temperature 
evolution of the VO 2 film. Notice that in these equations 
the only unknown parameter is the Kapitza constant ax- 
All the other quantities are known from direct measure¬ 
ments. P(Tc) is obtained from the profile of R{T) across 
the MIT, independently from the dynamics of R after 
the photo-induced insulator-to-metal transition. Also, 
the relation between the specific heat L and T is fixed by 
general and fundamental results^®^^^ . 

While these equations can in general be solved only 
numerically, some qualitative understanding of the de¬ 
cay time T can be gained if we make some simplifications. 
Let P{Tc) = 1/(v^o'tJ exp{-(Tc - {Tc)Y/\2a'f^^)}. At 
a temperature T the insulating volume fraction is given 
by ? 7 /(r) = i [1 - erf ((T- (Tc))/(v^crTj)] ■ Then as¬ 
suming that the pump pulse is strong enough to drive 
the entire film into a fully metallic state at t = 0, the 
time-dependence of the insulating volume fraction can 
be approximated by a simple exponential form 7?/(t) = 
1 — Ae~^l'^. In this case, an expression for the tempera¬ 
ture can be obtained through the relationship r]j{T{t)) = 
rji{t). Furthermore, assuming that the substrate temper¬ 
ature Tg does not change with time, and the latent heat 
L to be temperature-independent, we can calculate the 
decay constant r: 


r 


Cd 


{CmPM + LpavP{To)) 
o'x Tq — Tg 


+ To J 


(16) 


where the constants (7 > 0, and tq can only be deter¬ 
mined by solving the full system of equations (12)-(15) 

It is interesting to note that despite its many limita¬ 
tions, Eq.(16) captures many important qualitative traits 
of the actual relaxation process. For example. Figure (9) 
shows a plot of the decay constant as a function of ar^ 
obtained solving the full system of equations (12)-(14) for 
two different values of the average critical temperature, 
and same initial temperature Tq = 360 K. It is easy to see 
that the decay time r follows the linear trend predicted 
by Eq. (16) in the limit (Tg — {Tf}) S> \/2aT^- Similarly, 
an exact solution shows the inverse dependence of t on 
the Kapitza constant, t cx a~f^, as shown in Figure 10. 

The relation (16) is another important result of our 
work, it shows how the characteristic time of the recovery 
dynamics is related to the properties of the VO 2 films. In 
particular it shows the novel result that r grows linearly 
with the standard deviation of P{Tc). ax^ can be 
reliably obtained from the profile of R{T) across the MIT. 



dr, (K) 

FIG. 9. Full numerical calculation of the dependence of metal 
state decay constant r on for two different values of the 
sample average critical temperature (T,), and Ts{L) — 280 K. 
The initial temperature Tq = 360 K is such that the sample 
is initially fully metallic, and (Tq — {Tf})/ {-\/2aTf) ~ 9. 


V. EFFECT OF INHOMOGENEITIES ON THE 
RELAXATION DYNAMICS OF THE 
PHOTO-INDUCED MIT 

Using the theoretical approach described in Sec. IV we 
can obtain the time evolution of the optical reflectivity 
R through the MIT, as well as explain the significant 
difference in relaxation timescales between the two VO 2 
samples considered. In all the numerical calculations we 
assume C/, p/, Cm, Pm to be equal to the bulk values 
for insulating and metallic VO 2 , see Table II. 

The initial film temperature is fixed by the pump fiu¬ 
ence taking into account the Gaussian profile of the pulse 
and the fact that some of the heat is lost by the film dur¬ 
ing the time interval [0, to] for which our analysis does 















FIG. 10. VO 2 /AI 2 O 3 metal state decay time t dependence on 
the Kapitza constant ak for (D) = 64.7 nm, ao = 38.5 nm, 
substrate temperature Ts{L) = 310 K, and fluence (f) = 
8 mJ/cm^. The red dots correspond to numerical calcula¬ 
tions, and the dashed line is given by r oc 


VO 2 heat capacity insulating phase 

0.656 J/(g K) 

heat capacity metallic phase Cm*‘'^ 

0.78 J/(g K) 

density insulating phase 

4.57 g/cm® 

density metallic phase pm'^^ 

4.65 g/cm® 

thermal conductivity insulating phase 

3.5 W/(m K) 

thermal conductivity metallic phase 

6 W/(m K) 

bulk latent heat i,(Bnik)47 

51.8 J/g 

Ti 02 heat capacity 

0.686 J/(g K) 

density ps'^® 

4.25 g/cm® 

thermal conductivity 

8 W/(m K) 

AI 2 O 3 heat capacity 

0.779 J/(g K) 

density p^®^ 

3.98 g/cm® 

thermal conductivity Ks®^ 

30 W/(m K) 

V 02 /Ti 02 absorption coefficient at 800 nm a^® 

0.01 nm"® 

VO 2 /AI 2 O 3 absorption coefficient at 800 nm a^® 

0.0076 nm"® 


TABLE II. Parameters of VO 2 and substrates. 


not apply, t = 0 is time at which the pump pulse hits 
the VO2 film and tg = 10 ns for V02/Ti02 films and 
to = 0.5 ns for VO2/AI2O3 films. 

As discussed in Sec. IV, uk is the only unknown pa¬ 
rameter. For the case of V02/Ti02 samples, by fitting 
one of the curves for the dynamics of the reflectivity, 
we find ax = 1100 W/(K cm^). We find that all ex¬ 
perimental curves are well approximated assuming the 
same value for the Kapitza constant, see Fig. 11 (a). For 
the case of VO2/AI2O3 the characteristic timescale of 
the recovery is much shorter than for V02/Ti02 sam¬ 
ples. As discussed in Sec. Ill the two samples have very 


different inhomogeneities: is almost 4 times larger 

in VO2/AI2O3 than V02/Tid2. All other things being 
equal, Eq. ( 16 ) implies, see Fig. 12 , that r should be 
larger in VO2/AI2O3 than in V02/Ti02, the opposite of 
what is observed experimentally. We are then led to con¬ 
clude that ax in VO2/AI2O3 must be much higher than 
in V02/Ti02. Figure 13 shows the measured evolution 
of R for the VO2/AI2O3 sample for a fixed value of the 
fluence (j) and substrate temperature, and the theoretical 
curves for this case that we obtain using the distribution 
P{Tc) obtained for VO2/AI2O3 and two different values 
of ax- We see that by choosing for ax the same value 
used for V02/Ti02, there is no agreement between the¬ 
ory and experiment. By setting ax = 13000 W/(K cm^) 
we obtain excellent agreement. Indeed, all the experi¬ 
mental curves R{t) for VO2/AI2O3 are well approximated 
by the theoretical results assuming ax = 13,000 W/(K 
cm^). 


a) 



FIG. 11. (a) Time evolution of reflectivity after the photo- 
induced MIT for V 02 /Ti 02 for three different Ts{L) and (p = 
9 mJ/cm^. The solid curves correspond to the theoretical 
results, and the dashed curves correspond to the experimental 
results. For the three theory curves we use ax = 1100 W/(K 
cm^). Panel (b) shows the corresponding insulating fraction 
time evolution. 

Figure 14 shows the time evolution of the VO2 film 
and substrate temperatures (close to the interface) for 
the VO2/AI2O3 film, panel (a), and for the V02/Ti02 
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odUd) 

FIG. 12. Dependence of the V 02 /Ti 02 metal state decay 
time constant r on ct_d for two values of (InD), as defined 
in Eq. (5), Kapitza constant <jk = 1100 W/(K cm^), sub¬ 
strate temperature Ts{L) = 280 K, and initial fluence (j) = 9 
mJ/cm^. 



tins) 


FIG. 13. VO 2 /AI 2 O 3 reflectivity time evolution after photo- 
induced MIT for (j> = 7.5 mj/cm^. The red dots correspond 
to the experimental result. The dotted curve correspond to 
the theory with ax = 1100 W/(K cm^), and the solid curve 
corresponds to ck = 13000 W/(K cm^). 


film, panel (b), using the parameter values summarized 
in Table II. It helps to qualitatively understand the differ¬ 
ences in the thermal relaxation between the two samples. 
Due to the lower values of the Kapitza constant, thermal 
energy stays more concentrated near the V 02 -Ti 02 in¬ 
terface, keeping the temperature of the VO 2 film above 
Tc longer. 

To investigate the temperature dependence of the 
thermal relaxation we repeated the measurements 
while changing the base substrate temperature of the 
V02/Ti02 sample. For these measurements the sample 




FIG. 14. Film and substrate temperature time evolution. For 
sapphire (a), Ts{L) = 310 K, and for rutile (b), Tb{L) = 
280 K. 


was placed inside a cryostat, and cooled down to temper¬ 
atures Tg{L) between 260 K and 298 K. The results of 
these measurements, along with the theoretical calcula¬ 
tions, are shown in Figure 15. We again observe a good 
semiquantitative agreement between theoretical and ex¬ 
perimental results. Also, note that the simple expression 
for the decay constant r Eq. (16) captures the overall 
decay rate drop at lower substrate temperatures Ts{L). 

We point out that all the theoretical curves are ob¬ 
tained using the fixed set of parameters shown in Table II. 
As mentioned above, the only unknown parameter that 
enters the theory is aK- In the results presented above 
aji was fixed to a single value for each film, and this 
value was then used to obtain the results for a range of 
experimental conditions with different substrate temper¬ 
atures and pump fluences. For example. Fig. 4 shows an 
excellent agreement between the experimental measure¬ 
ments and theoretical calculations across the entire range 
of pump fluences, limited on the lower end by our ability 
to reliably detect the variation in the probe reflectivity, 
and on the upper end by the damage threshold of our 
sample (pump fluence > 40mJ/cm^). 



FIG. 15. Dependence of metal state decay constant r on 
fluence and substrate temperature for V 02 /Ti 02 . 
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VI. CONCLUSIONS 

In conclusion, we have presented a combined experi¬ 
mental and theoretical study of the long timescale recov¬ 
ery dynamics of VO 2 films following an ultrafast photo- 
induced insulator-to-metal transition. We have devel¬ 
oped a theoretical approach that is able to properly take 
into account the effect of inhomogeneities. The main 
results of our work are: (i) The derivation of the general¬ 
ized heat equation (9) that properly takes into account 
that during the recovery, due to the inhomogeneities, 
only fraction of the sample is undergoing the metal-to- 
insulator transition and correctly tracks the evolution in 
time of the metallic (insulating) phase; (ii) The clarifica¬ 
tion of the connection between the temperature depen¬ 
dent profile (i?(T)) of the reflectivity across the thermally 
induced MIT and its dynamics after a photo-induced 
insulator-to-metal transition; (iii) The approximate re¬ 
lation, Eq. (16), between the characteristic time of the 
recovery dynamics and the parameters of the film, in 
particular to the standard deviation of the distribution 
of critical temperatures as extracted from R{T); (iv) The 
ability of our theory to describe, using a fixed value of the 
Kapitza constant, the recovery dynamics for different val¬ 
ues of the substrate temperature and pump fluence. By 
changing the pump fluence the characteristic time of the 
recovery can be changed, experimentally, by two orders 
of magnitude: our theory is able to account for such a 
change. 

The theoretical approach that we present is general 
and can be used to describe the dynamics (in the adia¬ 
batic limit) of inhomogeneous systems across a first order 
phase transition. The approximate relation between the 


characteristic time r and the parameters of the system 
shows that r is directly proportional to the width of the 
thermally activated transition. This result allows to es¬ 
timate the recovery time of VO 2 films solely on the basis 
of a measurement of R{T) across the MIT. 

Assuming that variations of the size of the grains form¬ 
ing the films are the main source of inhomogeneities, us¬ 
ing very general and fundamental relations between the 
grain size and the grain’s critical temperature, we have 
been able to obtain the distribution of the grain sizes. 
In particular, we have been able to estimate the aver¬ 
age grain’s size and its standard deviation. We find that 
the calculated average grain’s size is in remarkable semi- 
quantitative agreement with the one obtained from XRD 
measurements. For systems in which inhomogeneities are 
mostly due to variations of the size D of the grains, our 
analysis provides a way to obtain the size distribution 
P{D) from the temperature dependent profile of the re¬ 
flectivity across the thermally induced MIT. This could 
be very useful considering that P{D) is a very challenging 
quantity to obtain via direct measurements. 
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